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ABSTRACT 


Based on a vertical motion model derived by Martin (1964), and 
applied by Jarrell (1967), a method is developed for producing diabatic 
peace function fields at 50 mb from independently computed fields of 
digthurte and adiabatic vertical motion. | Vertical motions over and above 
the adiabatically cero” pent are found to be of major significance 
in the diabatic heating function calculations. The resulting heating 
function fields are shown to be gompat te with existing synoptic trends 
at 50 mb and above, for the dat Stine B-5 April 1963. A proposed 


circulation Linkage involving ozone-heating areas is advanced, which 


fits all of the mean zonal computations of diabatic heating and of 


vertical motion into a consistent pattern. 
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L. Introduction 


Atmospheric heat budget analysis has become increasingly 
prominent in the environmental sciences literature of the past decade. 
The availability of lower stratospheric data has permitted the objective 
ye-examination of previous heating effect theories which were, unavoid- 
ably, restricted to global and hemispheric generalizations - beyond the 
local effects. As the quantity and reliability of such data increases, 
concepts will gradually lend practical perspective to synoptic scale 
variations in the development of numerical-prognosis models involving 
diabatic effects. 

Some of the notable investigations on adiabatic heating in the 
‘stratosphere during recent years are those of Budyko (2), Davis (4), 
and Ohring (11). Examples of more specialized investigations leading 
to synoptic scale applications are those of Jarrell (7), Martin (8), 
Paulin (13), and Pressman (14). The present study further focuses 
the attention os effects of lower stratospheric heating and vertical 
motion phenomena by investigating the 50 mb diabatic heating function 
for a time near the spring reversal during April 1963, 

Background for the present study began in 1964 when Martin 
developed a method for obtaining 12-hour mean velocity potential 
fields. This method was applied to 300, 200, and 100 mb levels for 
the period 1-5 April 1963 in order to study a time near the spring 
reversal. 

The data utilized by Martin became available from NMC in late 
1963, and did not extend above 100 mb. The significance of Martin's 


diagnostic technique for the present study lies in its application 


£3 





toward defining vertical motion fields in terms of geopotential, 
which are generally compatible with the quasi-geostrophic modeling 
equations [see Arakawa (1) ] . 

Working with Martin, Jarrell produced 50 mb vertical motion 
fields using a method which will be discussed below. These omega 
fields are based directly on Martin's ‘&-fields for 1-5 April 1963, 
and constitute one of the basic inputs to this study. 

Recently, 30 and 50 mb fields of both temperature and geo- 
potential (official analyses of the Air Weather Service) have become 
available for the same time and in the same format as the data used 
by Martin; that is, the octagonal 1977-point NMC grid. This new data 
will be utilized, in addition to Jarrell's omega fields, to compute 
diabatic heating function fields, which can then be used to determine 


zonally averaged heat transfer mechanisms. é 


2. Vertical motions at 50 mb 
Martin's process involved the numerical solution of the balance 


equation 


a 2 a, xe > 
ve(eon) +2] 5 Se Gy) | FP 


to produce nondivergent, rotational stream function fields in terms 


of geopotential, (0) = gz. By iterative solution of 


ee 


gat AN (nl) 
ge Bh = |e ETM - ME” 





15 





Martin solved for 12-hour mean potential fields for O600GMT and 
1800GMI at 100, 200, and 300 mb levels for 1-5 April 1963. The 
term in the brackets of (2) was the forcing function, and was known 
at each map time from the solution of (1). 

In an effort to determine vertical motions at 50 mb, Jarrell 


obtained smoothed divergence fields from 


VV HVE © 


computing the right side of (3) from Martin's stored X-fields. In 
(2) and (3), the superior wavy bar stands for the time centered value 
of the indicated parameter. The Laplacian operator in (3) was applied 


to the nine point grid. shown in Fig. Ll, with the usual format 


Ve v at mn Us +X% eae Th 4Xo (4) 


in terms of the nine point differencing mesh of Fig. 1. In (4), m 
is the map factor of the polar stereographic projection, Next, the 


continuity equation 


en 
Vp WV = dp (5) 
was integrated to compute vertical velocity fields at 300, 200, 100, 
and 50 mb. In order to do this, it was necessary to fit ow /2p by 
a Lagrangian interpolating quartic polynomial, It was further con- 


sidered appropriate to include the assumption ow/op = 0 at p=0, 


The resultant W  -fields were then expressible in terms of linear 








combinations of ie at 100, 200, and 300 mb, and maximum values 
at 50 mb as high as 10 cm/sec resulted. Concerning the validity of the 
assumption 3w/ap = 0, there is considerable evidence, particularly 
by Pantske (12), that both horizontal wind components during the 
stratospheric reversal are decreasing and in process of reversing sign 
on a hemispheric scale. Since the reversal normally occurs in mid- 
April, the assumption ow/ap= 0 seems appropriate to this general 
time, as noted by Jarrell. 

Therefore, Jarrell's work has produced thermodynamically un- 
restricted GW-fields at 50 mb which are suitable for comparison to 
independently computed Uyp~fields, This comparison will then lead 


to a solution for diabatic, or (Wp-fields. 


3. Diabatic W-equation 


The thermodynamic equation may be written in the form 


i ee AL, Be ane 
oe Ee ae Sie (6) 


Upon expansion of the time derivative and use of the geostrophic 


approximation for advection, one then obtains 





ei oT ‘ gt RT 
Cp = 20 aNg VT + w (2 CpP. | 7) 


By employing the hydrostatic consistency condition, 
T = -r ad {dnp ) for temperature in terms of geopotential, 
and using a potential temperature equivalent for the term in 


parentheses, eq. 7 becomes 


L7 





_peauog og RTp de 
ag | aby ae SPS)! 


This form has been chosen in order to utilize a static~stability 
parameter introduced by Holl (6) which he defines in the form 

Ss Skee 

H © op (9) 

The use of On is felt to be appropriate because it has been found 

by Fleet Numerical Weather Facility (FNWF), Monterey, California, to 

provide vertical consistency between standard analysis levels, while 


permitting horizontal consistency with randomly introduced wind data. 


This static-stability parameter may also be written é - 





~ | 
at RT _| FO ad K=B ao) 
="RP Op Cp ~ | dnp p = 


Since there is a second derivative in (10), Oy must be treated as a 

constant over a depth of three successive pressure levels in finite 

differencing about p = 50 mb. With the finite differencing pressure 

interval specified as ln P3/P, = ln 30/50, the third pressure level 

is taken as p = 83.3333 mb, with Py = 50 mb, and b= 30 mb, Since 

mean temperature, and the thickness equivalent (to be denoted X) 

in (15), below, will be shown to be proportional to one another; 

the choice of the 83.3 mb level is required to "match" 30 mb for an 

equal logarithmic interval across 50 mb.) Geopotentials at 83.3 mb - 


will be obtained by an interpolation process to be discussed later 


in this section. 











For the purpose of expressing Ox in terms of geopotential, 
g= g(r) is selected in the form of a quadratic polynomial in 


terms of In p, as follows 
P 2 
d = O,,+ adn 56 +b(Ln&) (1) 


By substituting p = 30 and p = 83.3 and adding and subtracting 


ek (7) 0 (® 12 
° 2D, ° >a) 4 


where 
Ve- b,- Dos vo i b,-26,,+ Ooas ie 


In (13), vo and vo are the first and second order finite 


difference operators, in the vertical, 





The thickness tendency equation. Because of the greater "noise" 


in the stratospheric temperature fields available to us, as compared 
to fields of geopotential, it is advantageous to use geopotential 
thickness (mean temperature} fields at 50 mb. This is a slightly 
different stability treatment than that of Craig et al (3) in the 
study of vertical motions associated with “explosive stratospheric 
warming". However, the foregoing development obviously calls for 
the geopotential thickness field approach. 


Eq. 8 may now be written 


o Se ae 
Stump 8 V3inp* B= ~ ee 


(14) 





Multiplying each term by p § (in p) and integrating from P, * 83.3 to 
P3 = 30 yields 
P 


pee, Sin sfaritefacer f POSdmp «5 


Eq. 11 suggests a definition of "thickness" parameter as follows 


Ps 
r= {(bpsie 
A 
so that B P, 
x =o SMnp = RT] P § Lop “6 
P3 P3 


The last form of the integral (17) has the value of 
RT(p, - P3) = 53,333 RT. Thus X = 53.3333 RT and is a measure of 
the mean temperature, and conversely of |the geopotential thickness. 


This confirms the statement made just above eq. 11. From eq. 11, 


3 @ fain p may be written 


ae = np 2, np Ln 54.2 ales 
= Qai\+2b (Ln $) (18) 


Then by (16), with Py » Py and P3 as previously specified, and with 


a and b as given by (12), X is written 





! = * bso 50 | 53, P Dey 
7 2 Ga 2 Whee 53,3 * aaa (solnt -833Ln%) 165 


= 34,7297 ®.. + 34.9478 $..- 69.6775 P¢3.3 


in units of gpm(mb), or ee) hb , For purposes of deriving a 
value for omega at 50 mb it is useful to consider K as "thickness" 
throughout this study. 

Now eq. 15 may be integrated to produce the diagnostic equation 
from which adiabatic-omega, GJ,» will result. Considering On (P) 
constant with respect to p in the layer Py to P3 (but not constant 


with respect to x and y), the result is 


é P3 
: = +450 ‘VX +6, ur, 53> ia a Q dp - 
P 


with X defined by (16). In (20), a superior wavy bar denotes a 


12-hour mean value. In finite difference form, (20) becomes 


MA 4+ Fad) (Bro)-hortes G05 =R(2) 53.3 (21) 


nd aad 
The solution for J and On is discussed below, 
Note that each term on the left hand side of (21) has units 
=] An. . =i 
of gpm(mb) (12-hrs) ~ so long as WW is in mb(12-hrs) ~. On the 
Lead 
right hand side, the column mean heating rate (Q) will be given units 


of j(12-hrs) (gm); so that RQ/Cp still has units of tem) *(12-hrey 


Therefore, the right hand side of (21) becomes (0.28704) (53.3) (Q/Gp) ’ 





j mb(gm) ~~ aCe ny ee But since g = 980.665 cm(sec) “, (from Ref. 18), 
ill 


it follows that one joule = 101.9716 gpm gm. Thus to obtain con- 
sistent units in both sides of (21), the right hand side must be 


multiplied by 101.9716, and eq. 21 becomes 


Ax +r F (G2o:X)~ suet G0 = |.5538xId Q (22) | 


in units of gpm(mb) (12-hrs)"!, 


The Jacobian term. The finite difference Jacobian of equations 


21 and 22 will have the usual form 


J ($, X) a (Sy OK q *X,) 7 (94- b) (Es X) (23) 


where the five point grid is formed from a subset of the nine point 
configuration shown in Fig. 1. The products formed by the J - 
operator have the units teow” mb. However, the term Bo is 


equal to 


mm 0.473 646 x 10 sec 
4 sing (mtr)* 





where m is, as before, the map scale factor! of the stereographic 
projection, Combination of the units of J in (23) with those of 


ar yhea- sate 


mb (4pm) sec _ q.¢oo.5 SS Salar 


(mtr)* 


Where meters and map factor appear together, meters will be denoted 
by mtr. 








fe 


2 
since 1 gpm = 9.80665 (mtr)/sec. The last statement merely redefines 
lL gpm in terms of MKS units. By transforming a to a 12-hour time 


interval, the Jacobian term has the full form 


2 f me 
OS5S01bS X lo (Fi), a (Zs. JX) (24) 


in units of gpm(mb) ‘ip nahi. 


Finally, because consecutive 12-hourly data is used, an average 


One 


Jacobian operator ( J ) is defined at each grid point. 


- At At 
(* =) +Ji (+ +4) 
z 


Ji 

= - (25) 
Ji t 

where At = 12-hours. 


The static-stabilit arameter. 





Eq. 22 contains a time averaged 


stability parameter ( Cn ) which may be expressed for each map time 


_yvd ve 
On Grape — 8 30a) Fh 


as 


consistent with equations 10 and 12; with K = R/Cp. Eq. 26 can then 


be reduced to the form 


G, = 40g, - 7460079, + 3.5527 Pte! ee 


which is the time averaged vertical stability parameter for the layer 


83.3 mb to 30 mb. The units of Gy are gpm. 
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Interpolation of g to the 83.3|mb level. Before the pre- 


viously developed expressions for Oy 


must be known. 


The synoptic levels of 


$ 83.3 


100, and 50 mb being 


and X can be computed, 


200, 


available; a parabolic fitting considered compatible with (11) is 


utilized. 


b (1n p) about p = 


6,,,(0=0 


100 mb, 


“Ca, 


Fioo 


Then, (28) can be expressed in finite 


to 


B54 = 0.2005 >. + 08/78 e 


Adiabatic Cp-solution. 


The resulting expression is 


+4(lo 


Equations 27, 


a Taylor expansion of 


$33) 


a) ee LY 


difference form and reduced 


(28) 


- 2.0480 d (29) 


19, 24, and 22 may 


now be assembled to produce the working equation for the solution of 


a : " 
pori (4.1120 $,,-7-6647 $_,t3: 


A. (34.7247 é, #40799, -69,67796 g, a) bs 


~2.f rity; 
ed X 
O.SOUWS KIO ( . J, (2xX 


or more simply 


_ AeX +0. 50165xI6" 


Wy 
mute 1 0uG On 


in units of suas) - 


Note that X 


5527 \ 00, = 
fs if 


(30) 


(31) 


neh j fi) Ji(eo% x) 


is used in two terms, that 


> 83:3 must be precomputed, and that 12-hourly averages apply 














throughout. The operator At is the 12-hour "time-differencer", 


and Z is the (dimensionless) contour height of the 50 mb surface, 


50 
It should be noted that the advecting wind is the geostrophic 
wind, rather than the quasi-geastrophic one. Also the static-stability 
(Gn) has been allowed to have space variability inherent in its 
definition, rather than held to a constant atmosphere value. Both of 
these energy constraint violations could assume major preportions in 

a prognostic procedure, but are minor in a diagnostic one. The use 

of Oy and of geostrophic advection were conveniently adaptable to 
the use of existing FNWF subroutines, 


On the octagonal boundary (U4) is taken to be zero at each 


boundary point. This is equivalent to choosing as boundary conditions 
ox = O and J} (2<0;3) =O 


at each boundary point. One grid interval within any boundary 
peint, both dX/2et and J Zso> X) are well defined, and a, is 


computable as at any other interior grid point. 


4, Heating function solution 


Eq. 20 may be written 


Gomes mr _ 100 R fi 


ated We 
For the adiabatic case, Q is zero, and G4 is zero. By using an 


unrestricted C& which is independint ty” computed, together with the 
ZNote the "time centered"Wis not defined as in (25), but results from 


(32). 
25 





just computed Gy , a solution for Wp is 


W- Wy = Wp 


As noted in Section 2, the value of W will be taken as that com- 
puted by Jarrell. By writing (32) for W and Gy , subtracting, = 


and using Ow from (27) we may write 


@ (2957) 12 = —hora(-2,)(4. a warg2s5i9,..) = 


A field of Gp is formed in units of mb (12-hrs)7* to compute the 
. be me a =r ail 
field of Q at 50 mb in j(gm) “(12-hrs) ~. 
_Recall that the right side of eq. 22 with qq in igey) “(tae ~ 
was converted to units consistent with those of (33) by the multi- 


plicative factor 1.5538 X 10°. Thus temporarily consider (33) as 


5538 X Io = 1.02164 Gp On (34) 





and obtain, upon simplification 


ore, —3 Pv and } 

Qi =-O6575XI0 Oy Wp (35) c 

ss 

. : 

In (35), Ox is the expression within the last pair of parentheses ; 


of (33), and, of course, is time averaged for each point of the grid, 
similar to the operation of (25). With (Cp also in mb(12-hrs), the 


notation Q, infers units of j/gm/(12-hrs), Finally, it should be 


noted that eq. 33 requires the assumption that the sum of the local 





change and advective change in geopotential thickness is the same 


whether (32) is written for W or for Wy, . In support of this 





assumption, a heat balance study for 1-5 April 1963 applied to the 
100 mb level by McClosky (10) indicated that vertical heat divergence 
was larger than heat transfer by (1.) horizontal advective and (2.) 
storage processes by one and three orders of magnitude, respectively. 
If this is also true at 50 mb, then (35) computes only the that 
is expressed by the vertical motion transfer, which appears to be the 


dominant contribution. 


5. Procedures and initial results 
The data grid and W-boundary effects. The FNWF 1604 (modified) 
CDC computer was used for all data processing, computations, and grid 
printing discussed herein. Figs. 2, 3, 4, 5, and 6 for 1 April 1963 
at O600GMI are presented as being typical of any of the ten map times 
processed in this study. Although the computer input data was initially 
in the form of a 47 X 51 point octagonal grid, the boundary of the 
computational octagon was restricted by four grid spaces to one of 
39 X 43 mesh size due to anticipated boundary value excursions in- 
troduced by (W (discussed by Jarrell). As a result, approximately 
10° of peripheral grid coverage was lost. However, the loss is not 
considered too significant for the purposes of this study, due to the 
sparse data area involved, and to the remaining grid area, which 
covers the Northern Hemisphere north of (approximately) 18° latitude. 
The above mentioned boundary excursions are evident in the 
GO -field, as seen in Fig. 2 (after Jarrell), and also in the 


aad ee tod 
fields of Wp and Q (Figs. 4 and 6, respectively) which depend 


directly upon Jarrell's 50 mb GO -fields. More recent investigation 
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PIG. 2. Twelve-hour average map of w vertically centered 
at 50 mb for O1 April 1963 at 0600cMT (mb/12-hrs). 
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FIG. 3. Twelve-hour average map of us vertically centered 
at 50 mb for 01 April 1963 at O600GMT (mb/12-hrs). 
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FIG. 4.  Twelve-hour average map of @Jp vertically centered ee 
at 50 mb for 01 April 1963 at 0600cGMT (mb/12-hrs). 








, ed 
FIG. 5.  Twelve-hour average map of Oy vertically centered 
at 50 mb for 01 April 1963 at O600GMIr (meters). 
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FIG. 6. Twelve-hour average map of Q vertically centered 
at 50 mb for 01 April 1963 at O600GMT (j/gm/12-hrs) 
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by Martin? (9) promises a marked decrease in the magnitude of this 
boundary effect on WW), but the revised data fields were not avail- 
able for this study. The boundary conditions for the , solution 
of this study simply set a= 0 at all boundary points and this bound- 
ary condition led to no apparent excursions in UW. However, the 
CO) -fields (see Fig. 3), computed by the method of Section 3, 

have also been "restricted" to the smaller grid of Fig. 5, within 
which Jarrell's 50 mb QW-values have been considered valid. It 
is to be recalled that Gy, is computed from eq. 31 within the 47 
by 51 octagon. 

Away from the boundary zone, the contoured C&)-dependent 
maps for 1 April 1963 at O600GMT have a regular appearance with a 
similarity in scale patterns. However, OO and cd, field centers 
do not coincide. Even with the anomalous (/-values near the 
boundary, a primary goal of this paper was that of determining the 
principal meridional heat transfer mechanisms by using, where possible, 
the zonally averaged quantities [Wp], [Ou] ,[Q] . In the averaging 
around latitude bands most of the spurious boundary influences are 
averaged out. 

Presentation of initial results. Fig. 4 shows merely the 
point-wise subtraction of the fields of 0, from w , GE Bde. 77, 
and the result is the "diabatic" WO, -field. Fig. 5 shows the 


ODA Oe 
field of Oy for the initial time in the sequence. Here Oy is 


3the proposed approach involves setting the boundary inflow-outflow 

component of Vy equal to the divergence of the geostrophic wind at 
all points centered one row in from the 47 X 51 octagon in a manner 
similar to that of suppressing “spurious anticyclogenesis" in the 
W -field. 
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FIG. 7. 


Five-day mean map of a vertically centered at 
50 mb for 1-5 April 1963 (mb/12-hrs). 
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FIG. 8. Five-day mean map of ain vertically centered at 
50 mb for 1-5 April 1963 (mb/12-hrs), 
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FIG. 9. Five-day mean map of ®Wp vertically centered at 
50 mb for 1-5 April 1963 (mb/12-hrs). 
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FIG. 10. Five-day mean map of ro vertically centered at 
50 mb for 1-5 April 1963 (meters). 








Be ih: 


- 
oN 
& 
\ 


FIG, 11. Five-day mean map of a) vertically centered at 
50 mb for 1-5 April 1963 (j/gm/12-hrs). 
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computed at each grid point by eq. 19, but only that portion of the 
On field within the truncated grid has been reproduced, as 
discussed above. Comparison of Fig. 5 with Fig. 18 [50 mb chart from 
the analysis series of the Free University of Berlin (15)] shows a 
remarkable coherency between locations of higher (lower) values of 
On and the location of warm-ridge (cold-trough) areas; which is a 
reasonable result. 

Fig. 6 shows the first field of computed i’ in j/gm/(12-hrs), 
centered vertically at 50 mb. The computation of Qi; at each grid 
point is based on eq. 35 with the values of ( Oa );, and ( 0p dij 
recalled from the taped records of the applicable data fields. 

Similarly, fields of co 4 OO, 5 Hy > and Sy were also stored 
at appropriate grid locations, for each of the other nine map times. 
Thus a five day mean map was generated for each parameter. These 
maps are labeled a , Wy » Wo, Gy and Q in Figs, 7-11 
respectively. The superior bar denotes a five day time-mean. In 
the following sections, the bracket notation enclosing any of the 


above fields will denote its zonal mean. 


6. Zonal means 

Computational procedures. Zonal averages by 5° latitude bands 
were computed for each parameter listed in Table 1, A FNWF subroutine 
was utilized for obtaining zonally averaged values centered at 20°N, 
25°N++-90°N. Then all ten fields of each parameter were averaged 
together to produce the results seen in Table 1. The "average" for 


the pole point is simply the computed average value for that point 


weighted with the four surrounding points. A separation of 0600GMT 
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is the diabatic vertical velocity, 
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Zonal means at 50 mb by five degree latitude bands for 1-5 April 1963. 
parameters are five day averages: [W@W] is the vertical velocity, [&4] 
the adiabatic vertical velocity, [Wp] 
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All 


[Gy] is the vertical stability parameter, and [Q] is the diabatic heating 


function. 





zonal means from those at 1800GMT to permit a possible diurnal 
investigation was not attempted since the sample was not considered 
large enough. 

Inspection of a preliminary version of Table 1 revealed single- 
point sign reversals of [W] ,(p], and [@] at 30°N and 50°N. 
Whether these sign reversals indicated actual transitory surges, or 
computational outgrowths of unrepresentative data, is not known. 
Since these were small perturbations superimposed on larger latitudinal- 
scale patterns, a one-pass spherical smoother (described below) was 
applied to the preliminary table. 

In order to smooth an arbitrary function of latitude, U =U (9), 


the Laplacian was applied to the one dimensional array ui), 


[Ui], = [TU], +4 ¥, [TI(aAgy (36) 


where Vs is the spherical finite difference form of the Laplacian. 
The azimuthal and radial terms of the full spherical Laplacian are 
not applicable, so that eq. 36 contains only a zonal differencing 


term of form 


=: ee = Cae 


a” a* cos 


Thus the values of Table 1 have been adjusted according to 


[Vi] ne [Ud + z (vel -2[Ui)+ [U%)- avi] {Ui-))tan 0: Fa (37) 


where “i increases toward the pole, and the smoothing angle between 
successive i values is S§1f/(180) radians. The "end point" values 


at ¢ = 20°, and at ¢ = 90° were not subjected to smoothing, although 
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FIG, 12. Graph of five-day mean stability parameter by five degree latitude bands for 
1-5 April 1963 (meters). Horizontal line at 1971 meters is the standard 
value, using eq. 27 and standard. geopotential heights from Ref. 18. 
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FIG, 13. Graph of five-day mean vertical motion parameters by five degree 
latitude bands for 1-5 April 1963 (mb/12-hrs). 
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FIG. 14, Graph of five-day mean diabatic heating function by five degree 
latitude bands for 1-5 April 1963 (j/gm/12-hrs). 


as with all other entries in Table 1, they were time averaged, 

In order to interpret the most essential aspects of Table 1, 
[On] ; [uy] and [Wo] » and (Q] are graphed on a latitude scale in 
Figs. 12, 13, and 14, respectively. Note that singular departure 
points are not present in Table 1, or in the associated graphs; and 


meridional patterns are evident in all of the tabulated parameters, 


7. Zonal mean eddy heat-flux and mean-cellular heat-flux convergences 
Since Q is expressible by means of eq. 35, it follows that the 


zonally averaged [Q] is expressible as 


[Q] = -K[Og vo] — K [oy* wo*] K= 0.0575 X160° “sf 


at each latitude, and for each map time, The first term on the right 
is the contribution of the meridional circulation flux convergence 
(denoted Q. }) common to the latitude band. ‘The second term is then 
ascribed to eddy flux convergence of heat (denoted Qe ), resulting 
mainly from synoptic scale systems. Since [Q] and -K [Su] Lo] = [Qe] 
were known for each map time and latitude, the eddy flux term [Qe] = -K 
[Oy*tu] in (38) was computed as a function of latitude and map time. 
Then the values of [Qe] and (Q¢] were time averaged and smoothed in 
the same manner as the Table 1 parameters. The results ciutedmy 
zonal mean values are presented in Table 2 under the headings; "cellular" 
and “eddy flux convergence contributions to [0] ". oe [o] values 
from Table 1 are reprinted in Table 2 to facilitate comparison with 
[@] and [Qe] . 
It would also be possible to perform a time mean analysis of 


[a] in terms of the time mean values of the right side of (38), and 
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-,08 +,03 +.19 +,.22 +.06 +.03 +,05 +.01 


-.09 +~-,03 +,.22 +,32 +.12 


-.05 +.04 +,21 +.25 +,06 +.02 +.06 +,02 


[Q. ~.08 -.08 -.03 -.01 


-.08 -,02 +,22 +.31 +.13 


-,02 -.03 -,00 +,01 ~,01 -.OL1 -.O1 -,O1 -.00 +,01 -,O01 


TABLE 2, 


Zonal means at 50 mb by five degree latitude bands for 1-5 April 1963, 


All 
parameters are five-day averages: [Q) is the diabatic heating function, 


[Q] is the cellular flux convergence contribution to [Q] , and [Q] is the 
eddy flux convergence contribution to [Q] . 
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FIG, 15. Graph of five-day mean diabatic heating function paraubbenp: (Qj is 


the cellular flux convergence contribution to [Q] , and{Q.] is the 
eddy flux convergence contribution to[Q] . 


of their transitory flux transports. However, the transitory aspect 
of these terms has been subordinated to that of determining the two 
time-averaged distributions of vertical flux convergences, which have 
been graphed as functions of latitude in Fig. 15. 

The results of the [Q.] computation indicate that it is 
generally negative in sign, with one minor exception [ near 50°N 
(Fig. 15) | . However, [Qe] is, for the most part, an order of 
magnitude smaller than [Qa] so that for diagnostic purposes we may 
equate [Q] and [Q] . Note that the horizontal flux convergence 
could be involved in the [Q] distribution, but the assumption was 
made that no substantial change occurs between the diabatic and 
adiabatic cases; an assumption made also by Paulin (13). Hence the 

[Q) distributions assessed here are those additive contributions 
arising from transformations in the vertical. However, for this 
particular period, as contrasted with cases of the onset of explosive 
warming, maximum values of the horizontal transport flux convergence 
terms have been found to be an order of magnitude smaller than those 


resulting from vertical exchanges [ see McClosky (10) | “ 


8. Interpretation of the results 

A standard comparison for atmospheric heating studies is the 
excellent research of Davis (4). Since Davis' computations were 
primarily based on radiation and climatological air mass data, a com- 
parison with the present study can only be made if its results are 
regarded as deviations from the mean April computations of Davis. 


Values presented by Davis for ten degree latitude bands from 


20°N to 70°N, are taken as mid-band averages, converted to j/gm/12-hrs, 





: DAVIS 


25° 35° 45° 55° 65 75 85 
- 068 -.179 - 280 - ,350 - 403 
THIS STUDY 
25° 35° 45° 55° 65° a 85° 
sed +.038 +.176 +041 +,001 -.005 +, 269 


TABLE 3. Averaged April 1963 heating rates at 50 mb in ten 
degree latitude bands (j/gm/l2-hrs). 
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and listed in Table 3 from 25°N to 65°N. These values are based on 

a column centered at 40 mb, but are felt to be sufficiently representa- 
tive of 50 mb for the purposes of this comparison. The values listed 
in Table 3 for "this study" are taken from Table 2, using a three 


point averaging technique of the form 


[ely — §o% BAC: +é cos ZAI; + COS Din | Qin (39) 


a COS Pj, +2.c0S Dj + cos Piz, 


where gi increases poleward in five degree increments. 

It is seen that the (30-day average) values taken from Davis 
call for cooling in all latitude bands, and with magnitudes steadily 
increasing poleward. In comparison, the (five day average) results 
of this study show a larger cooling rate in the subtropics, and warming 
throughout most of the remaining hemisphere. The subtropical cooling 
rate, however, is representative in sign only, due to the & boundary 
value effect at this latitude. The relatively small cooling rate 
centered near 70°N is quite definite (see Figs. 14 and 15), but is 
terminated poleward by the strong warming centered near 85°N. 

Further investigation reveals that an extrapolation of Davis’ 
infrared cooling and solar heating rates apparently undergo a reversal 
by mid-April, presumably near the 35 km level [see Figs. 2 and 3, 
and Table 4; Davis (4) | . However, it would be impossible to infer 
from such an extrapolation, the latitudes of preferred energy input 
at levels below 35 km, since the climatological smoothing eliminates 
any possible consideration of radiation-dynamics in the stratosphere. 


A hemispheric radiational study to explain the 50 mb warming effects 
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Deviations of the April 1963 absolute topography from 
the 50 mb mean for 1949-53 (geopotential dekameters). 
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for the period of this study (Figs. 14 and 15) was too prodigious a 
task to perform with the same degree of detail utilized in the thermo- 
dynamic computations employed herein. Moreover, many data gaps exist 
on a daily basis, which cannot be smoothed over as in a purely climat- 
ological approach, 

Therefore, the alternative of a judicious synoptic comparison 
has been attempted, in order to find out how much the data sample of 
1-5 April 1963 differs from synoptic stratosphere climatology, and 
whether these differences are compatible with the 50 mb (Q] distribu- 
tion resulting from this study. 

The synoptic charts to be utilized are those from the series 
published by the Free University of Berlin (15), (16), (17). Fig. 16 
illustrates the mean 50 mb height anomalies for the month of April 
1963. The developing Arctic anticyclone was more prominent than 
normal, as was the subtropical belt of high pressure. However, a 
considerably greater anomaly is obvious in the strong negative height 
deviations which dominate the Central Asian and North American areas. 
These negative anomalies are particularly characteristic of actual 
conditions during the period of this study, as evidenced by Fig. 18 
and 19. The southward movement of the subpolar vortices should be 
noted, as well as the general broadening of the associated trough 
areas. These trough areas have major southwest and southeast exten- 
sions from the corresponding coastal areas of both continents in the 
vicinity of 45° latitude. This feature is evident on the mean April 
map of Fig. 17, but is particularly prominent during the data period 
of 1-5 April. For reasons which will be discussed below, a major 


significance of these trough areas is their vertical extent above 50 mb. 
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FIG. 17. Mean absolute topography of the 50 mb surface 
for April 1963 (geopotential dekameters). 
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FIG. 18. 50 mb height analysis for 01 April 1963 at 0000GMT * 
(geopotential dekameters). 
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FIG. 19. 50 mb height analysis for 06 April 1963 
at OOO0GMI (geopotential dekameters), 
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FIG. 21. 10 mb height analysis for 01 April 1963 at 
OOOOGMT (geopotential dekameters). 





$7 


Figs. 20 and 21 may be considered representative of the period 1-5 
April at 30 mb and 10 mb respectively. The cyclonic circulation 
patterns referred to at 50 mb are even more pronounced at 30 mb ’ 
and 10 mb. 

Therefore, it is apparent that the previously noted areas of 
negative height anomalies are also areas of cyclonic baroclinic activity, 
both of which have considerable vertical extent above 50 mb. The pre- 
dominence of these conditions during the time of this study has 
important implications regarding the heating function distribution of 
Fig. 14. Also, it may be noted that this persistence of winter regime 
stratospheric conditions is during a month in which the transition to 
the summer regime easterlies is normally in process. It should be 
noted that the spring reversal is known to have been considerably 
delayed in 1963. 

The development of Section 3 has relagated horizontal transfer 
processes to a secondary role. Energy generation by diabatic influ- 
ences will be considered the primary mechanism to be investigated J | 
relative to an explanation of the heating function distribution | 

s 
determined by this study. It is apparent from Figs. 13 and 14 that 
(4) and [ud] have an opposite sign relationship. Therefore, energy 
source (sink) regions may generally be located by negative (positive) 
values of [wp] . Fig. 13 shows major regions of [a] > 0 near 45°N 
and 85°N, whereas the region of the subpolar vortices is associated 
with [aj}<o . The effect of high level ozone irradiation in polar 


latitudes is considered to play a significant role in the inter- 


pretation of the 50 mb diabatic heating distribution resulting from 
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Schematic model of proposed vertical circulation patterns for 1-5 April 1963. 
"SH," is the subtropical high, "I.E." means trough extensions, "S.V." means 


subpolar vortices, and "A,H," is the Artic high. Vertical motions vectors 
represent the relative intensity of &# from Fig. 13 (&# =U +We ). 





this study. Fig. 22 schematically presents a proposed vertical 
circulation model with which the various synoptic and dynamic con- 
siderations discussed above may be combined. ’ 

In Fig. 22, proposed locations of ozone-radiative heating 
maxima are at A, and Ay. The primary maximum at A, is well documen- 
ted by the ozone climatology of Gétz (5). The secondary maximum at 
Ay is based on the high level advection of ozone into southwest and 
southeast extensions of large-scale, broad-based troughs near 45°N. 
This synoptic feature is evident in Fig. 21, and is presumed to exist 
at 40-45 km also. Therefore, distinctive cyclonic circulation 
patterns are transporting ozone away from the Gétz high level maximum 
to the region of 45°N, Radiative warming then produces in situ 


expansion at A,, though probably of smaller intensity than at A 


2 Ta 


The forced outdrafts at A, and Ay maintain the proposed thermal - 
cellular circulations of Fig. 22. These patterns can be shown to 
be consistent with all the vertical motion and diabatic heating 
results of this study. Notice, in particular, that the 50 mb warming 
( faq> O) at 45°N indicates the trough extensions to be filling 
features. This is substantiated by Fig. 23. 

It must be understood that Fig. 22 is speculative in the region 
above 50 mb, and cannot be substantiated by this study. However, 
the linkage with high level ozone advection, and with existing 
synoptic circulation patterns, as well as complete agreement with 
the results of this study, shows an intriguing consistency which may 
warrant further investigation, Nevertheless, Fig. 22 must be regarded 
only as a hypothetical circulation pattern which is consistent with 


the results produced by this study at the level of 50 mb. 
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FIG, 23% Mean absolute topography of the 50 mb surface 
for May 1963 (geopotential dekameters), 
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9. Concluding remarks 

From the results that have been presented, the relationship 
between vertical motion and diabatic heating is best demonstrated 
by &p . Recall that Ud, was used as a convenience for producing 

U4» . Furthermore, it should be noted that both Wp and Wy vary 
widely from &/ , but in a manner which suggests that the importance 
of Wr», is linked with the importance of external heating effects. 

In polar latitudes where irradiation heating at 50 mb was minor, Us 
was of much larger magnitude than G/p) . This interpretation is 
consistent with the hemispheric importance of Wp as indicated by 

this study, and with the proposed cellular circulation pattern of 

Fig. 22. Previous studies have produced similar Gp values as 

found here, but have evaluated the diabatic effect on adiabatically 
deduced vertical motion fields to be relatively small [see Paulin (13). 
However, by reason of time and sample differences, and differences 

in the development of Gp» values, it is felt that the present results 
are not contradictory. 

The Martin QY-field method has led to 50 mb U/-fields and 
from there to Q-fields which describe synoptic scale diabatic heating 
effects in the lower stratosphere. No vertical velocity data is 
known with which to directly compare the W-fields produced by 
Jarrell, However, the Us fields of this study afford an order of 
magnitude check which, it is felt, generally support Jarrell's 
findings. Furthermore, an averaging technique similar to (39), when 


applied to the [Ww] values of Table 2 yields a hemispheric average 


of (essentially) zero. But most significantly, these M-fields 





have circulation features which can be reasonably interpreted, and 

are consistent with synoptic maps for the same time period, More- 

over, the Q fields resulting from this study show similar consistencies, 
These results are considered sufficiently promising to justify 

further development of the method which produced them, and application 


of this method at still higher levels; for example, at 10 mb. 
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